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Abstract: According to Dijkgraaf and Vafa the effective glueball superpotential of 
the M = 1 supersymmetric QCD coupled with an adjoint chiral multiplet is given 
by the planar amplitude in the j? expansion of a matrix model. It was shown that, 
when the M = 1 supersymmetric QCD is coupled with fundamental chiral multiplets 
as well, the effective glueball superpotential is modified by the disc amplitude of the 
generalized matrix model. The diagramatic computation of this disc amplitude is 
fairly involved for the multi-cut solution. Instead we compute it with recourse to the 
complex analysis of the hyperelliptic curve. The result is given in series of the gluino 
condensation Si. The explicit computation for the generic multi-cut solution is done 
up to order S 3 . It is systematic so that it can be extended to higher orders. 
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1. Introduction 

Dijkgraaf and Vafa discovered a close relation between the M = 1 supersymmetric 
QCD and a matrix model. They claimed that the glueball superpotential of the M = 
1 supersymmetric QCD coupled with an adjoint chiral multiplet can be computed by 
the -k expansion of an associated matrix model]!]]. It was confirmed by subsequent 
works by many people. The effective theory of the M = 1 supersymmetric QCD 
given by such a glueball superptential is called the Dijkgraaf- Vafa (DV) theory. In 
[0] they developed a remarkable technique to compute the planar diagram amplitude 
of the matrix model. Namely the computation, which was normally done by the 
perturabtive sum of the planar diagrams, was equivalently replaced by the analysis 
of the hyperelliptic curve on the Riemann surface. This revealed important aspects 
underlying in the M = 1 supersymmetric QCD like the duality 0, the integrable 
hierarchy, etc, and attracted a renewed attention 0, |4| to the Seiberg-Witten theory 
for the N = 2 supersymmetric QCD[|5|. 

In H |7|, [8], H the relation between the the Af = 1 supersymmetric QCD and a 
matrix model was generalized for the case in which the Af = 1 supersymmetric QCD 
is coupled with fundamental chiral multiplets as well. The presence of fundamental 
chiral multiplets modifies the glueball superpotential of the Af = 1 supersymmetric 
QCD. It was shown that the modification can be captured by computing the planar 
diagrams with boundary (discs) of the matrix model which is generalized accordingly. 
The disc amplitude is the perturbative sum of such disc diagrams. It was computed 
to all orders of the gluino condensation S for the 1-cut solution0, but merely to 
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order S 2 for the n-cut solution[||. The computation to higher orders for the n-cut 
solution is extremely involved. In this paper we will discuss that the disc amplitude 
for the n-cut solution may be computed to higher orders in a rather systematic way, 
if one has recourse to the complex anlysis of the hyperelliptic curve. Indeed the 
planar amplitude for the 2-cut solution was calculated along this line to order S 3 in 
PHI . This calculation to order S 3 was generalized for the n-cut solution in [|11|] later 
on. Our aim of this paper is to generalize the calculation of the planar amplitude in 
[LlJ to that of the disc amplitude. 

The paper is organized as follows. We give a short review on the DV theory 
coupled with fundamental chiral multiplets in section 2. The disc amplitude of 
the associated matrix model is given in terms of the (hyper)elliptic function. As a 
warming-up it is evaluated for the 1-cut solution in section 3. The evaluation is exact 
for this case. For the multi-cut solution it can be done only in series of the gluino 
condensation S. In section 4 the disc amplitude for the 2-cut solution is evaluated 
to order S 3 . In section 5 this evaluation is generalized for the n-cut solution. 



2. A Short Review on the DV theory 

We consider the Af = 1 supersymmetric U (N) gauge theory with the superpotential 



Nt 



w($, Q 7 , Qi) = trv{§) + J^i^QiQ 1 + Qi^Q 1 ), 



(2.1) 



i=i 



Here $ is a adoint chiral multiplet and Q 1 (Qi) are Nf chiral multiplets in the fun- 
damental representation. V(<&) is a polynomial scalar potential of degree Mo given 
by 



V'(x) — (x — a\){x — a 2 ) • • • (x — a Mo )- 



(2.2) 



When the U(N) gauge group symmetry is broken to UiU(Ni) with Y^=i = this 
gauge theory is described by the effective action 



W, 



eff 



^l + F 1 (S) + 2 7 rir j2Si 



(2.3) 



i=i 



for the gluino condensation Si. It was claimed that F2(S) and Fi(S) are the planar 
and disc amplitudes in the expansion of the partition function of the associated 
matrix model 



Z = I dcjxiq 1 dqi exp 



1 



W^q^qj) 

9s 



exp[-J>;*F x (S)' 

X<2 
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Here 4> is a M x M matrix and q 1 ^ 1 ) are M dimensional vectors. W is the super- 
potential given by fl2.1|) so that the matrix model is gauge- invariant under U(M). 
The gluino condensation Si is identified with the 'tHooft coupling g s Mi, when the 
gauge symmetry U(M) is broken to Hill (Mi) with Y17=i = ^ ■ ^he planar and 
disc amplitudes, respectively given by F 2 and F 1; are calculated by an alternative 
method. Namely we think of the bosonic matrix model 

f r 1 

Z= dcftexp trV(4>) 

J L 9s 

Then we obtain the saddle point equation for the eigenvalues distributed with density 
0(A) 

y = V'(x) +2 [ dX-^- = 0. (2.4) 
J x — X 

Here y is the quantum deformed force of the matrix model. For the n-cut solution 
it may be written as a hyperelliptic function such that 



y=(x- nx)(x - Hi) ■ ■ ■ (x - fi Mo ~n)V(x ~ Ai)0 - A 2 ) • • • (x - X 2n )- (2.5) 

Comparing the polynomial part in the expansion of ( J2.4Q and ( |2.5| ) at x = oo gives 
the zero-points a a , a = 1, 2, ■ • • , M as functions of /i's and A's. We consider the DV 
differential ydx defined with such a hyperelliptic function. For the n-cut solution the 
gluino condensation is given by 

Si = — / ydx. (2.6) 

The planar amplitude [0J was given by 

<9F 2 



55i 

while the disc amplitude M, O by 



A 

ydx, (2.7) 



Af/ A 

F ' = \il\l ydx-V(A)+V(-m I )}, (2 



7=1 J 



with a cut-off parameter A. 

We want to calculate this disc amplitude for the n-cut solution and compare it 
with the result obtained by the expansion. For simplicity we shall do it simply 
assuming that M = n. 
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3. The 1-cut Solution 



As a warming-up we start with the simplest example, i.e., the matrix model with 
the scalar potential 

V{x) = 

which is coupled with one chiral matter multiplet. We consider the DV differential 

ydx = ^V'(x) 2 - A^Sdx, 

on the Riemann surface with one branch-cut. Note that the branch points are chosen 
so that ( |2.6| ) is satisfied, i.e., 

if 

2m I 



4S 



The disc amplitude ( |2.8|) is exactly calculated as 

n , m S 1 . r- — - _ „, 1 + \J\ - AaS ni m 

F 1 (S) = ---—(VT=toS-l) + Slog + Slog 

with a = This is in complete agreement with the result obtained by the 

expansion 0. The log-term which is linear in S is a non-perturbative part and may 
be absorbed by the last term in the effective action ( |2.3|) . A similar calculation for 



the 1-cut solution was done for a variant of the matrix model in [O . 



4. The 2-cut Solution 

Next we study the matrix model in which the scalar potential is given by 

V\x) = (x-a l )(x-a 2 ). (4.1) 

In H they calculated Feynmann diagrams dipicted in Fig. 1 to find the disc amplitude 
to order S 2 : 

Nf 2 r S 2 1 2 S S- 1 S 2 i 
Fi(S) = - [ K ^J- i — - 2 g Xlah;,/,, + 2 A^jf 2 ! ■ (42) 

Here Aj(a) = Uj?n a ijj fu — m i + a ii a ij — a i ~ a j- in this section we shall evaluate 
the formula Q2.8p for the 2-cut solution to order S" 3 and show that it indeed reproduces 
the result (14. 21) to order S 2 . To this end we consider the DV differential 



ydx 



y((^-7i) 2 -p 2 )((^-72) 2 -^)^, (4.3) 
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Figure 1: The Feynman diagrams to order S 2 corresponding to the three terms in (|4.2|), 
The solid double line and the single dashed line are propagators for <fi and q 1 respectively. 
The solid and dashed double line is a propagator for ghosts. 



on the Riemann surface with two branch-cuts [71 — a±, 71 + «i] and [72 — ct 2 , 72 + a 2 ]- 
Here y may be written in the form Q2.4| ). Hence comparing the polynomial part in 
the expansion of ( |2.4|) and ( fOl ) at x = 00 yields the relations cti(p 2 , 7): 

a 1 + a 2 = 7i+72, (4.4) 
(«! + a 2 ) 2 + 2«ia 2 = 7 2 + 7 2 2 + 4 7l72 - p 2 - p 2 . (4.5) 

With the DV differential ( |4.3|) we will calculate the disc amplitude ( |2.8| ). A suitable 
technique to this end was developed for the calculation of the planar amplitude in 
p"T|j . The calculation goes following the steps: 

1. evaluation of the period integrals 

1 pi+Pi 

Si = — ydx, (4.6) 

which provides explicit formulae Si(p 2 ,j) for gluino condensation. 

2. inversion of the formula Si(p 2 ,7) which provides p 2 (S,~f) in series of Si with 
coefficients made out of 7^. 

3. inversion of the formulae a.i(p 2 ,j), given by ( [4.4[ ) and ( [4.5| ), which provides 
7i(p 2 , a) in series of p 2 with coefficients made out of «j. 

4. evaluation of the DV differential ( |4.3|) in series of p 2 with coefficients made out 
of if 

5. expansion of p 2 and 7^ in the above result in series of Si with the coefficients 
made out of a, by iterating the inversion formulae obtained at Step 2 and 3. 



Step 1 



After the shift x' = x — 71 we expand the DV differential ( ^4.3|) in series of p\ as 



ydx = c m pl m (x' + mY'^Jx' 2 - p 2 dx', (4.7) 



m=0 
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where 712 = 71 — 72 and c r , 



r(m-|) 
mir(-i)- 



Write (x' + 7i2 



,l-2m 



for m > 1 as (2m — 2)-ple 



derivative of (x' + 712) 1 by 712. Then it takes the form 



ydx 



P 2 r 



m> 



^(2m-2)!^ K d ll2 



9 ^2m-2 V /x ' 2 ~ P? 

^' + 7i2 J 



dx'. (4.8) 



We integrate the DV differential of this form along the Ai-cycle surrounding the 
branch-cuts [71 — p%, 71 + pi]. We find 



pi 
-pi 



c/x'(x' + 712) yp? - x' 2 



J7T 



-7l2Pi, 



(4.9) 



for the first piece and 



,2 _ x /2 



Pi 



+ 712 



CT (-\/7i 2 2 - Pi +712) 



(4.10) 



for the second piece. Putting them together in (fO|) then we get the expression for 
the gluino condensation Si 



q 1 2 PlP2 PlPj(Pl+Pl) , n(n 8 X 

^1 = T7l2Pi - "J + <->(P )• 

4 8712 327j 2 



(4.11) 



Similarly we get for the gluino condensation S 2 



S 2 = -7712P2 + v 

4 8712 



P?P2 , P?P2(Pl+Pl) 



+ 



327i 3 2 



+ 0(p 8 



(4.12) 



Step 2 



Inverting the relations ( [4.11 ) and ( 4.12|) we obtain 



Pi 



45i 8^5 2 245 1( S 2 (^i - 5 2 



712 



7i 4 2 



7l2 



+ 0(S 4 



2 _ 45*2 8SiS 2 245 l 1 ff 2 ('S l i — 6*2) n / Q 4 
P2 — 2 7 *~ 



712 



7i 4 2 



7l2 



(4.13) 
(4.14) 



Step 3 



Inverting the relations ( |4.4| ) and ( |4.5| ) we obtain 



7i 



72 



-(«! + a 2 ) + i^( ai -a 2 ) 2 -2(p? + p|), 

Pi + P2 (pi + pI? 



2a 



12 



4a 3 



0(p fe 



12 



-(«! + a 2 ) - 2 V ~ tt2 ) 2 ~ + ^2) 
(ft + 4±£f + W r H«! + 0(A 



2a 



12 



4a? 2 



(4.15) 



(4.16) 
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Step 4 



We expand the DV differential ( |4.8|) in series of pi as well 



ydx 



V'(x) 



112, PI 



p\ 



1 _hfi__ 4 



2 x' x' + 712' 8 \r 



^ + 7i2) s 



P2 1 



•}+ 



x' 3 (x' + 712) 3 4712 x' x' + 7i 2 ' 



Jo 



16 x' 4 (x ; + 712) 4 J 16 x' 5 (x' + 712) 5 

4^2 



} 



+ 



+ 



1 ,P\PI 



{ 



PlP2 1 ^rPlP2 



P?P2 



I6712 x 13 (x' + 712) 3 167^ x ' 2 ( x ' + 712) 2 



} 



167i 3 2 



V x' + 712 



+ 0(p 8 ) 



dx' . 



(4.17) 



Here V'(x) is given by 



V'(x) = x' 2 + ll2 x > --(pl + pl), 



which is nothing but the potential ( |4.1|) owing to ( |4.4|) and ( |4.5|) . After shifting the 
variable x' back to x we write (|4.17| ) exressing pf and 73 in terms Si and a« with 
recourse to ( f4.13|) ^( }4.16|) . In particular the terms of order p 2 and p 4 in (|4.17|) read 



712, Pi 



— \ — — — ) — 2{S\ 0-S1S2 p-S\S 2 (Si — S 2 )} 

2 x' x' + 712 a'{ 2 cn" 2 x — a\ 

9 3fl 1 

+ 2{S 2 + —S.S2 + —S 1 S 2 (S 1 - S 2 )} 

af 2 «i2 x — a 2 



4{4-Si(Si - S 2 ) + -irS^Sl - 9S X S 2 + AS 2 2 )}- ' 



a\ 2 a\ 2 -■- (x - ai) 2 

+ 4{^ r S 2 ( 1 S 1 - S 2 ) + ^-S 2 (4S 2 - 9S 1 S 2 + 3S 2 2 )} . 1 

«i2 «i2 {x-a 2 y 

8 5i(S'i — S'2) 2 8 S 2 (Si — S 2 ) 2 

I" ^ 7 To ~T 



8 l x' 2 



Pi 



x' + 712) 2 



ot\ 2 (x — «i) 3 a\ 2 (x — a 2 ) 3 
} = 2{^l + ^S 1 (2S 2 1 -3S l S 2 )} ' 



«12 "12 



X — Oil) 2 



+ 2{4" + -4^(3^2 - 2S 2 2 )}- 

a\ 2 a\ 2 (x — a 2 ) 2 

8 S 2 (St-S 2 ) 8 SziSt-Si) 



+ 



712 r Pj_ P\ 

C L <3 / / 



a 12 ( X ~ a l) 3 a 12 ( X ~ a 2) 3 

- 2 4 _ 1 



S 

x' 3 (x> + 7i 2 ) 3_} = + ^fe^ 1 ~ 2S2)} \x - «i) 3 

-2{^l + ±S 2 2 (2S 1 -S 2 )} 1 



"Q!l2 «12 



x - a 2 f 
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p\p\ , 1 _ 1 ■ 

4712 V x' + 7i 2 ' 



12 - ff 2 ) 12 ^(fr - S 2 ) 

a 12 ( X ~ a l) 4 a 12 ( X ~ a 2) 4 ' 



1 



a 



12 



a 



12 



H ^-SiS 2 (Si — S 2 ){ 



+ 



x — a.\ x — a 2 
1 



a 



12 



(x — ai) 2 (x — a 2 ) 



With these the disc amplitude Q2.8Q is evaluated as 



TP I Q\ ST [ C 1 m / + "1 C 1 m / + "2 

Fi(5) = - 2^ log 1 S 2 log ■ 



1=1 



A 



A 



+ (4-( 5 i 2 - 2S1S2) + 4-(2S? - 9SfS 2 + 65 1( S 2 2, 
Va 2 2 a\ 2 

+ (^r(Sl - 25i5 2 ) - 4-( 2 ^2 3 - 9^1^ + QSfS 2 ] 



1 



12 



12 



m_f + ai 
1 

m/ + a 2 



+ 



S ' + -^-(2S 3 - + 2S 1 S, 



2«i 2 a\ 2 



Sj + ^-(2Sl - + 2S 2 S 2 



2ai2 «f 2 



(mi + at) 2 
1 

(m/ + a 2 ) 



2 r 1 1 

+ 3*[ 2 V 2S! ~ 3S2lS2) (mi + ai )3 " (2 ^ 2 ~ fo7 + 

^1 



(4.18) 



+ 



+ 



mj + a 2 ) 3 . 



St 



2a\ 2 I (mj + «i) 4 (m/ + a 2 ) 



0(S 4 



The log-terms linear in S is a non-perturbative part which may be absorbed by 
the last term in the effective action (|2.3|) . The terms of order S 2 are in complete 
agreement with the result ( |4.2|) by the 4 expansion. 



5. The n-cut Solution 



We shall generalize the calculation for the 2-cut solution to that of the n-cut solution. 
That is, taking the matrix model in which the scalar potential is given by 



V'(x) = (x — azi)(x — a 2 ) ■ ■ ■ (x — a n ). 
we consider the DV differential 



(5.1) 



n 

ydx = Y[[(x- 7i) 2 - ptjdx, 
\ i=i 



(5.2) 



on the Riemann surface with n branch-cuts. Again we calculate the disc amplitude 
following Step 1 to Step 5. 
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Step 1 and Step 2 



The inversion formula pf(S, 7) was given in the expansion form[l l 



Pi 



1 - 



j 



A^V 2A,( 7 ) j£ Jij7ik ' -^^(7) 



(5.3) 



in which 



A i(7) = Y[lij, 



Jij 7t 7j • 



As it will be clear, the expansion up to order S 2 is sufficient for the calculation of 
the disc amplitude to order S 3 . 
Step 3 

We expand y, given by fl2.4p and Q5.2| ), at x = 00. Comparing the polynomial 
part of both expansions yieldsjl I 



2 Z) «i = Z) [(7i + A)' + (7<- Pi 



i=l 



i=l 



1 < / < n. 



(5.4) 



We shall express 7« in terms of p and a by solving this set of equations recursively. 
First of all expand fl5.4| ) in series of p 2 as 



z °i = z N + 9 7i- 2 P? + ; i\-' P \ + o(/) 



i=l i=l 

Assume that 7$ is further expanded in the form 

7, = cu + °? ] + ^ + 0(p 6 



(5.5) 



(5.6) 



with a; the term of order p n . Putting this in (5.5) we get the recursion relations 

r (™) 
for 07 



0. 



1 

zK v f ) +^ / - i )«^ 



(5.7) 



z 

1=1 



+ ^(/-l)(/-2)(Z-3)a|- 4 A 4 ' 



0, 



and so on. Note that the identity 



AAa) ' 



1 j ' 



(5- 



(5.9) 
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for the quantity 



V'(x) 



Y[(x-ai). 



(5.10) 



Using the expansion 



H) becomes 



1=1 



M'atf- 1 , 



By means of this identity the recursion relations (|5.7| ) and Q5.8|) can be solved for 



(n) 

07 as 



cr 



(2) 



(7 



(1) 



1 71 

VpH'fa V 



(5.11) 



2 Ada 



3=1 



24A,(a 



1 n 

VF [,, '"f«-)o 4 



(5.12) 



3=1 



Here pW'(x), pW"(x), PW"' (x) are derivatives of the quantity ( p,10| ) with respect to 
x. Putting ( |5.11 ) and ( 5.12j ) in ( |5.6| ) we obtain the formula generalizing ( |4.15 ) and 
( f4.16| ) for the n-cut solution. 
Step 4 

Expanding y given by the DV differential ( |5.2|) in series of p 2 we get 

.2 1 .4 ! p 6 



2/ 



^L2x-7i 8(a;-7i) 3 1Q (x - %) 5 



n 



p 2 ip] 



+ 



1 



(s-7,)(x-7 i ) 16(x-7 i )(x-7 i ) 



P [i] (x) 

P [i %) 



(5.13) 



(x - %)(x - Jj)(x - -f k ) 



P m (x) + 0(p 8 ), 



with 



P®(x) = Y[{x-Ti), 

n 

P m (x) = H (z-y). 
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When P®(x),pM(x) and P^ k \x) are expanded at x — 7i,7j,7fc in Taylor series, it 
takes the form 



y = V\x) + J2J2 



(*-7«)" 



+ 0(p% 



with 



aii(p 2 ,7) 



Pi 



P [i] ""(%)-^P [1] "(%)-^P [ H%) 



384 



16 K ' 2 

A n 2 





167?- 4 7ij J 






G2i(p 2 ,7) = 


6 4 "42 

96 u ; 8 V ' 1 ^ I67?, V ' 


G3i(p 2 ,7) = 


6 4 n A 

^ W (7*)+Ei6 ipW1(7i) 


G4i(p 2 ,7) = 


-fp ] '(7*), 


a5i(p 2 ,7) = 


-^ [ %)- 



n ^4^2 



(5.14) 



The polynomial part V'(aj) was determined by the fact that y given by (|5.13 ) may 
be written in the form (|2.4|) as well. 
Step 5 



Finally we use the inversion formulae for p\ and 7$, given by (|5.3|) and (|5[ 
respectively, to find y in the form to order S 3 



5 n 



I/=l 1=1 



(5.15) 



with 



b u (S) = -2S h 



"2S^ — AS{Sj . Sf QSj Sj -\- SS^S*; Sf 3Sj Sj -\- 3SiS^ 



n 



L Aj(a)on 



+ 4- 



-4- 



A^A^af. 



Ai(a) 2 a|.a ifc 



A i (a)A i (Q;)Q;f i a; ifc 
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^^3SiSj ASfSj -\- 4:SfS k ASiSjSk ^S^S 2 + S^S^ 2SiSjS k 



A i (a)A i (a)af i a jfc 



A i (a)A j (a)a ij Q !j 2 fc 



+8 



-4 



C2 Q I C2 C 



2SiSjS k 



+ 



MS) 



hi{S) 

hi{S) 



n 



S{Sj ISiSjSk 
A i (a)A fe (a)a ii a| A , 

;S? + lOSfSi - AS 2 S k - lOS?Si + 



A i (a)A j (a)a ij a ik a jk 



A i (a) 2 a ij aikaii 



Y$>iSj — ASiSjSk + 12SiSjSi — ASiS^Si 



A i (a)A j (a)a ij a jk a j i 



-4 



C2c C2c 
<->i £>k — <->,; 



+ 



Q Q O O O O 



Ai(a)Ak(at)aija ik a k i Ai(a)A k (a)aija jk a k i 
2S? A r45 i (5 i - - 25,) 4S?(Si - 2Sj 



+ 0(5 4 



Ail a) 



A,(a)V 



n 

n 



ASi{S^-Sj)^i-2Sk 
Ai(a) 2 a i:j ai k 



*S!(Sj-s k ) 

Ai^Aj^aijCXjk 



+ 0(S 4 ), 



E 



ASf(2Si-3Sj 



AAaYa 



+ 0(S 4 



A, ; fa) 



Here the coefficient bu(S) is the exact result obtained by comparing the term tending 
to - at x = oo in (|2.4j ) and (|5.13[) . We would like to remark that the inversion formula 
(|5.3|) for p 2 to order S 2 is sufficient to evaluate the other coefficients b ui (Sys to order 
S 3 . Having explicitly known y in the expansion form ( 5.15Q it is easy to evaluate the 
disc amplitude Q2.8|) for the n-cut solution. We then see that for n = 2 it exactly 
reproduces the disc amplitude given by ( |4 . 1 8|) . Moreover we note that the evaluation 
of Q2.8|) to order S 2 gives the disc amplitude (|4.2|) in the extended form to the generic 
n, which was obtained by the expansion in ||. 



6. Conclusions 

To the author's knowledge the disc amplitude for the multi-cut solution was calcu- 
lated only the diagramatic expansion in the literature. The result was merely re- 
ported to order S 2 for the generic n-cut solution^]. The calculation to higher orders 
can be hardly done, because the diagramatic expansion becomes complicate as the 
degree of the polynomial scalar potential increases. To confront with this situation 
we exploited the alternative method, i.e., the complex analysis of the (hyper)elliptic 
curve. In section 4 the method was worked out to order S 3 for the disc amplitude 
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for the 2-cut solution. Dealing with the elliptic curve it is relatively easy to extend 
the calculation to higher orders. In section 5 the argument was generalized for the 
generic n-cut solution. By using a more systematic algorithm than in section 4, the 
disc amplitude for the n-cut solution was calculated up to order S 3 . There we have 
given the inversion formula for 73 to order p 6 . We would like remark that it was suf- 
ficient to give the formula to order p 4 for the calculation of the planar amplitude to 
order S 13 ^]. It was checked that the disc amplitude thus obtained indeed reproduces 
the one for the 2-cut solution obtained by the -k expansion in ||. The algorithm in 
this paper allows us to generalize the calculation of the disc amplitude for the n-cut 
solution to higher orders. 



Note Added: 

After posting the paper on the arXiv, the author was told that the disc amplitude 
for the n-cut solution was calculated to order S 3 by a different method in [|T3|| . But 



the result (4.16) in [13| is erroneous. After the correction it reads 



^-. x =i \^ T 1 f ^ Si(2Si — 3Sj) 2Sf 5SfSj + 2SiS, 

f=l i=l AK i e if j^i K6 U i e *j e if K i e ij e if 

QS 3 - 2QSfSj + lASjSj 2S 3 - 8SfS, + 2S t S] } " " 

+ ^ + ]+ hh 

by extracting the single sum contribution Ylij^A' ' '} fr° m 12k^i anc ^ ^j+i 

Ylk^j- Then this agrees with the disc amplitude given by (|2.8| ) with ( |5.15| ) in our 
paper, up to the double sum contribution X^=i Ylij^A' ' '}■ To check this, note that 



with the notation in 13 



n q N-l q N-l n n 

which follows from 

El dx / Si x dx 

, J, (x-e^W'ix) = ' tT^l (x-e^x-e^W'ix) = ' 

It is more involved to check the consistency between both results for the triple sum 
contribution Yli=i Tlij+i ^2k+j,ii' ' '}■ The above calculation uses the resolvent and 
the algorism developed in [|, |TJ|. It is simple, but allows us to evaluate merely 
counter integrals of the DV differential. On the other hand the calculation in our 



paper uses the other algorism ||11|| . It allows us to argue on the DV differential itself, 
expressing the branch points of the Riemann surface as functions of Si's. The latter 
algorism is important to study the Whitham hierarchy of the DV theory. 



13 
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